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ABSTRACT 

Using the IRAM 30-m telescope and the 15-m JCMT, we explore the value 
of para-formaldehyde (p-H 2 CO) as a tracer of density and temperature of the 
molecular gas in external galaxies. The target of our observations are the lobes 
of the molecular ring around the center of the nearby prototypical starburst 
galaxy M82. It is shown that P-H2CO provides one of the rare direct molecular 
thermometers. Reproducing the measured line intensities with a large velocity 
gradient (LVG) model, we find densities of n m ~ 7 x 10 3 cm~ 3 and kinetic 
temperatures of T^ n ~ 200 K. The derived kinetic temperature is significantly 
higher than the dust temperature or the temperature deduced from ammonia 
(NH3) lines, but our results agree well with the properties of the high-excitation 
component seen in CO. We also present the serendipitous discovery of the 4 2 — > 3i 
line of methanol (CH 3 OH) in the northeastern lobe, which shows — unlike CO 
and H2CO — significantly different line intensities in the two lobes. 

Subject headings: Galaxies: individual (M 82) — galaxies: ISM — galaxies: star- 
burst — radio lines: ISM — submillimeter — ISM: molecules 



1 Based on observations carried out with the IRAM 30- meter telescope. IRAM is supported by 
INSU/CNRS (France), MPG (Germany) and IGN (Spain). 



-2- 



INTRODUCTION 



Molecular gas is regarded as the fuel for star formation. Increasing evide nce for non- 



Paumard et al 



2006 



Klessen et al. 2007), and warm 



Rigopoulou et al.l 120021 ; iMauersberger et al.l 120031 ) 



standard initial mass functions (e.g. 
molecular gas in starburst galaxies (e.g. 
suggest that its physical properties may influence the star formation rate and the properties 
of the next generation of stars. Unfortunately, the physical properties of the molecular gas 
in external galaxies, in particular the kinetic temperature, are often not well constrained. 
The easily thermalized and optically thick CO J = 1 — > and 2 — >• 1 transitions could 
constitute a good temperature tracer, but the filling factor of extragalactic clouds is poorly 
constrained. Other commonly observed molecules like HCN and HCO + are good density 
tracers, but require an a priori knowledge of the kinetic temperature. The inversion lines of 
the symmetric top molecule ammonia (NH3) are frequently used as the galactic "standard 
cloud thermometer". However, in the disk of the Milky Way, th e fractional abundance of 
NH3 varies between 10 ~ 5 in hot cores (jMauersberger et al. 1987 ) and 10~ 8 in dark clouds 
( IBenson fc Myersl 119831 ). Thus, ammonia may trace preferentially a specific component of 
the molecular gas and the assumption of an approximately constant fractional abundance 
on linear scales of a few 100 pc is likely not valid. Other symmetric or slightly asymmetric 
top molecules may therefore be more favorable for extragalactic line studies. 



1.1. Formaldehyde as a Molecular Gas Tracer 



In this paper, we investigate the diagnostic properties of para-formaldehyde (p-H 2 CO) 
lines in extragalactic sources. Formaldehyde is formed on the surface of dust grains by 
successive hydrogenation of CO (IWatanabe fc Kouchill2002l ). released into the gas phase by 
shocks or UV heating and subsequently destroyed by gas-p hase processes. Variat ions in 
fractional abundance rarely exceed one order of magnitude. iJohnstone et al.l (120031 ) found 
only little variation in a variety of galactic environments ranging from cool protostellar 
candidates to hot photon-dominated regions and infrared sources with the latter showing 
the highest discrepancy, by a factor of 5 to 10. These higher abundances are likely caused by 
the release of H 2 CO molecules from grain surfaces in regions of massive star formation. As 
another example for the level of variability in the abundance of formaldehyde versus that of 
ammonia in active environments, the H2CO abundance is the same in the hot core and the 
compact ridge of Orion A, whe reas the NH3 abundan ce in the core is two orders of magnitude 
higher than that in the ridge (jCharnley et al.lll992h . 



Para-formaldehyde, which is a subspecies of the slightly asymmetric top molecule H 2 CO, 
has a rich mm- and submm-spectrum of transitions. Since the relative populations of the 
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K a ladders are generally governed by collisions, line ratios involving different K a ladders 
(inter-ladder ratios) are generally good tracers of the kinetic temperature. With the temper- 
ature known, line ratios within a single K a ladder, the so-called intra-ladder ratios, whose 
populations are mainly determi ned by collisional excitation and radiative de-excitation, sen- 
sitively probe the gas density (IMangum fc Woottenl Il993l ). The ^-doublet transitions of 
ortho-H 2 CO in the cm-range, lio — > In at A = 6 cm and 2n — > 2 12 at A = 2 cm, have been 



observed in a number of external galaxies since the 1970s (e.g. iGardner fc Whiteoak!ll974 



Seaquist fc Belli Il990l ) . H 2 CO lines in the mm-range we re detected in several galaxies in 



eludi ng the LMC, M 82 and NGC 253 since the 1990s (e.g. Baan et al.lll990l : Ijohansson et al. 



19941 ). but to date there are only three published studies of the excitation conditions of 



extragalactic molecular gas that utilize multi- line H2CO observations in combination with a 



non-LTE treatme nt of the radiative transfer (jrliittemeister et al. 1997 ; Heikkila et al. 1999 



Wang et all 120041 ) . 



1.2. Target Selection: M82 



The target for our exploratory investigation is the nearby prototypical starburst galaxy 
M82. As the most prominent galaxy in the northern hemisphere w ith a nuclear star- 
burst, M82 has been the target of numerous molecular line studies (e.g. iFuente et al.l 12006 



Seaquist et al.l 120061 ) . Most of the excitation studies focus on the central lkpc disk, where 
the molecular gas is concentrated in a circumnuclea r ring around the center of the starburst 



e.g. IGarcia-Burillo et al.ll2002l ; ISeaauist et a 



at smaller radii (Fig. [fl see also 



Mao et al 



112006!) with the more highly excited lines found 
"(hood ) for an overview of ring diameters). Note 
that the two lobes of the molecular ring in this highly inclined galaxy can be easily separated 
with single-dish telescopes and the observed molecular lines in the lobes have typical widths 
of only 100 to 150 kms -1 (e.g. ISeaquist et al.lll998l ; IGarcia-Burillo et al.ll2002l ). thus greatly 



(Mao et al. 


2000; 


Ward et al. 


2003) 



ggest that the broad data base of CO and 13 CO lines 
observed towards the molecular lobes of M 82 can not be fitted by a single- component large 
velocity gradient (LVG) model and that the majority of the molecular gas may actually be in 
a high-excitation state. The H 2 CO(li o — » In) transi t ion a t A = 6cm was detected towards 
the nucleus of M82 in absorption by I Graham et al.l (119781 ). whereas the H 2 CO(2n — > 2i 2 ) 



line at A = 2 cm was observed in emission, which was interpreted a s evidence for very dense 



(10 6 cm 3 ) molecular gas in the molecular ring (IBaan et al.lll9901 ). The latter study also 



reported the first detection of an H 2 CO emission line in the m m-range in M82. A first 



multi-line H 2 CO excitation study towards M82 was performed by iHiittemeister et al.l (11997 



see also Sect. II. ip using an LVG model that included 16 ortho-H 2 CO levels of the K a = 1 
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ladder. 



1.3. Line Selection 

Having developed a comprehensive LVG code for ortho-H 2 CO (40 levels) and para-H 2 CO 
(41 levels, see Sect. 14.11) . we searched for the set of transition lines best suited for extragalactic 
observations while trying to minimize observational uncertainties. The usually small filling 
factor of extragalactic observations and the decreased sensitivity of ground-based receivers at 
high frequencies limits the choice to the strongest lines, usually corresponding to transitions 
between low excitation levels. Among the H 2 CO lines that can be observed with present-day 
receivers, the para-H 2 CO transitions 3 3 — > 2 02 , 3 22 — > 2 2 i and 3 2 i — > 2 20 stand out by being 
close in frequency and at the same time strong enough for extragalactic observations. With 
rest frequencies at 218.22 GHz, 218.48 GHz and 218.76 GHz, all three lines can be observed 
simultaneously, if the receiver and spectrometer both offer a bandwidth of 1 GHz. That way, 
the observed inter-ladder line ratio H 2 CO(3 3 — > 2 02 )/H 2 CO(3 2 i — > 2 20 ), which is a good 
temperature tracer, will be free from uncertainties related to pointing accuracy, calibration 
issues or different beam widths. Due to the large line widths of extragalactic observations 
(> 100 kins" 1 ), the H 2 CO(3 22 — ► 2 21 ) line may be blended with methanol emission (see 
Sect. El) and thus be of limited use for the LVG analysis. At the very high densities (10 6 cm -3 



Baan et al.lll990l ) and kinetic temperatures of 50 K to 100 K inferred for H 2 CO in M82, our 
LVG code suggests a line ratio H 2 CO(3 2 i — > 2 20 )/H 2 CO(5 2 4 — > 4 23 ) of the order of unity, 
which sensitively traces the gas density. In addition, observing the H 2 CO(3 2 i — > 2 20 ) line 
with the IRAM 30-m telescope and the para-H 2 CO(5 2 4 — > 4 23 ) line with the 15-m dish of the 
JCMT results in data with nearly the same beam width (11" vs. 13") and thus minimizes the 
uncertainty inherent in line ratios derived from observations consisting of single pointings. 
Using a spectrometer with a bandwidth of at least 900 MHz, the H 2 CO(5 2 4 — * 4 23 ) tran- 
sition at 363.95 GHz can be observed simultaneously with the weaker and heavily blended 
para-H 2 CO transitions 5 42 — > 4 41 and 5 4 i — » 4 40 at 364.10 GHz as well as with the blended 
ortho-H 2 CO transitions 5 33 — > 4 32 and 5 32 — ► 4 31 at 364.28 GHz and 364.29 GHz, respectively. 
The non-detection of the H 2 CO(5 24 — > 4 23 ) line (see Section [3]) hinted at a much lower gas 
density regime than previously assumed and necessitated the addition of the relatively strong 
para-H 2 CO(2o 2 — > loi) line at 145.60 GHz to the selected transitions in order to be able to 
derive the density-sensitive intra-ladder line ratio H 2 CO(2 02 — > loi)/H 2 CO(3 03 — > 2 02 ). 



Here we present the results of our para-H 2 CO observations. To our knowledge, this is 
the first dedicated search for para-H 2 CO transition lines of the K a = 2 ladder outside of 
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the Galactic neighborhood (Milky Way and Magellanic Clouds). The distribution of the 
H 2 CO emission and its correlation to other tracers will be discussed in a forthcoming paper 
based on high-resolution (~ 4") ortho-H 2 CO data recently obtained with the VLA. After a 
description of the observations and the data reduction strategy in Sect. 2, we present the 
observational results in Sect. 3. In Sect. 4, we describe our LVG code and the parameter 
space covered, before we discuss our results and compare them to other molecular excitation 
studies. Finally, we report the serendipitous detection of methanol emission in one of the 
lobes (Sect. 5). Our conclusions are summarized in Sect. 6. 

2. OBSERVATIONS AND DATA REDUCTION 

For this study, we selected the para-H 2 CO transitions described in Sect. 11.31 to be ob- 
served at nearly the same spatial resolution with the IRAM 30-m telescope and the 15-m 
JCMT. The parameters of the observations are summarized in Table [TJ 

2.1. Observations with the IRAM 30-m Telescope 

For the observations at 218 GHz, the IRAM 30-m telescope with the heterodyne receiver 
array HERA, consisting of nine dual-polarization receivers arranged in form of a center-filled 
square with a pixel separation of 24", was pointed towards the southwestern lobe (SW lobe) 
of M82 at a 20 oo = 09 h 55 m 49 f 54, 5 200 o = 69°40'43'.'l during four nights in March/April 2005. 
The derotator optical assembly was set to keep the pixel pattern stationary in equatorial 
coordinates at an angle of 20° clockwise relative to the axis of right ascension, so that the 
central row of pixels was aligned with the major axis of the molecular ring with the two 
molecular lobes being covered by two adjacent pixels (Fig. [1]) and the third pixel pointing 
at a position 24" southwest of the SW lobe. In this arrangement, the remaining pixels 
of the array pointed towards six positions 24" northwest and southeast of the major axis 
(Fig. [2]). The spectra were obtained in a wobbler switching mode with a beam throw of 240" 
in azimuth. As backend, we used the WILMA autocorrelator, which provides a bandwidth of 
1024 MHz (512 channels) in each of its 18 units. Having tuned the receivers to a rest frequency 
of 218.48 GHz at a velocity of 200 km s -1 , each bandpass thus covered simultaneously the 
para-H 2 CO(3o3 — > 2 02 ), (3 22 — > 2 2 i) and (3 2 i — > 2 20 ) transitions (see Tabled]). Good winter 
weather conditions resulted in system temperatures of 200 to 400 K (T£) and a pointing 
accuracy of 3" or better. The chopper wheel calibration was checked by observing the 
C 18 0(2 — > 1) line at 219.56 GHz toward the reference source CLR2688 and the relative 
stability of each tuning was determined by observing the H 2 CO emission of NGC 2264 IR and 
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DR21 OH. Both checks suggest a calibration uncertainty of less than 15%. The temperature 
scale was converted to main-beam brightness temperature (T m b) adopting a main beam 
efficiency of B c r = 0.55 and a forward efficiency of F e ^ = 0.90. 

The para-H2CO(2o2 — > loi) spectra at 145.60 GHz were obtained at the IRAM 30- 
m telescope during two days in August 2006. Using the wobbler switching mode (200" 
throw) of the new control system (NCS), we observed each molecular lobe in M82 with 
receivers C150 and D150 tuned to the same frequency and each connected to two parts of 
the 1 MHz filterbank, resulting in a bandwidth of 512 MHz with a channel spacing of 1 MHz 
for each receiver. During the observations, the system temperature was 200 to 350 K and 
the pointing was accurate within 3" to 4". Daily observations of the galactic sources W51E1 
and/or W30H suggest a calibration uncertainty of 10% or less. For the conversion of the 
antenna temperature to main-beam brightness temperature (T m b) a main beam efficiency of 
B c s = 0.69 and a forward efficiency of F c s = 0.93 were adopted. 



2.2. Observations with the JCMT 

During four shifts in December 2004 and January 2005, the SW lobe of M 82 was ob- 
served with the dual-mixer receiver B3 at the JCMt| tuned to a rest frequency of 364.1 GHz 
and a velocity of 180 km s -1 . The observations were carried out in beam switching mode 
with a throw of 210". The DAS autocorrelator provided a bandwidth of 920 MHz with a 
spectral channel separation of 1.25 MHz in each of its two receiver channels, sufficient to 
cover simultaneously the para-H 2 CO transitions (5 2 4 — > 4 23 ), (5<u — > 4 40 ) and (5 42 — > 4 41 ) as 
well as the ortho-H 2 CO transitions (5 33 — > 4 32 ) and (5 32 — > 4 31 ) (see Tabled]). The system 
temperatures during the four shifts varied from ~ 750 K to ~1650K and the uncertainty 
of the chopper wheel calibration was about 15% as indicated by observations of the stan- 
dard source IRC+10216. The temperature scale was converted to main-beam brightness 
temperature using a beam efficiency of _B e g = 0.63. 



2 Program ID M04BC10; the James Clerk Maxwell Telescope is operated by The Joint Astronomy Cen- 
tre on behalf of the Science and Technology Facilities Council of the United Kingdom, the Netherlands 
Organisation for Scientific Research, and the National Research Council of Canada. 
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2.3. Data Reduction 



All data were reduced, converted to the main-beam temperature scale and analyzed 
using the GILDAS software package CLASS. Each individual spectrum was inspected for bad 
channels, standing waves or other baseline problems before a linear baseline was subtracted. 
The spectra at each position were then weighted according to their noise level, combined and 
smoothed to a common velocity resolution of lOkms -1 . The data obtained at 218 GHz with 
the newly commissioned WILMA backend required special attention, because a large fraction 
of the individual spectra showed baseline problems ranging from prominent DC offsets to 
subtle ripples in certain sections of the bandpass. As a first step, the spectra were sorted into 
different groups according to their baseline quality. In order to guard against more subtle 
baseline problems compromising the determination of line strengths, we then fitted the lines 
in each individual spectrum using the multi-line fit of the corresponding combined spectrum 
(Section [3]) as a template. When we selected only the spectra with the best baseline quality, 
representing ~ 40% of the data, the x 2 -tests yielded an agreement between the distribution 
of the fitted intensities and the expected Gaussian distribution at a 95% confidence level for 
each line in the spectra. The uncertainty of each fitted line intensity shown in Table [2] is the 
uncertainty in the mean of the Gaussian distribution er</> = aj/y/N with 07 the standard 
deviation of the Gaussian distribution and N the number of selected spectra. 



3. RESULTS 

We have detected the H 2 CO(2 2 — > loi) line at 146 GHz as well as the transitions 
H 2 CO(3 03 -> 2 02 ), H 2 CO(3 22 -> 2 21 ) and H 2 CO(3 21 -> 2 20 ) at 218 GHz in both molecular 
lobes of M82. This constitutes the first detection of H 2 CO(3 22 — > 2 2 i) emission from an 
external galaxy. We did not detect any H 2 CO lines at 364 GHz, probably because of the 
limited sensitivity of our observations. The spectra of the molecular lobes and the fitted 
Gaussian lines are displayed in Fig. [3j The basic assumption that the observed lines are all 
emitted from the same volume of molecular gas with the same average physical properties 
implies that all lines should have the same velocity profile, parametrized in a Gaussian fit by 
the central velocity and the line width. Therefore we fixed the velocity and line width of all 
spectral features to the values derived from simultaneously fitting all the lines in the 218 GHz 
spectra. The derived velocities and widths in the two lobes are in very good agreement with 
the results of a fit to the strong H 2 CO(3o3 — > 2 02 ) line alone. The parameters of the Gaussian 
fits and the resulting integrated intensities are summarized in Table [2j 

The 218 GHz spectrum of the SW lobe is reproduced well by the three H 2 CO lines near 
this frequency. The central velocity is 132kms _1 (local standard of rest, hereafter LSR); the 
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line width is 111 kms -1 . In the spectrum of the NE lobe, the central velocity is 301 kms" 1 
(LSR) and the line width is 131 kms -1 . The apparent shift of the H 2 CO(322 — ► 2 21 ) line 
relative to the central velocity indicates the presence of at least one additional emission line, 
which we identify as the methanol CH 3 OH(4 2 — > 3±E) transition at 218.44 GHz. The ratio 
of the fitted intensities of the highly blended lines H 2 CO(322 — * 2 2 i) and CH 3 OH(4 2 — > 3i E) 
can vary strongly with slight changes in the central velocity or the line width, which suggests 
that the fitting procedure is strongly influenced by the noise in the spectra. In addition, 
there might be faint HC3N(24 — > 23) emission at 218.33 GHz in the spectrum of the NE 
lobe. However, the effect of such an additional component on the other fitted Gaussians 
is well inside the observational errors and can be neglected. The integrated line intensities 
resulting from a simultaneous fit of all five emission lines are given in parentheses in Table [2j 
A simultaneous fit of all five lines to the spectrum of the S W lobe gives a stringent upper 
limit to the HCaN(24 — > 23) emission at this position. 

In the SW lobe, the two H 2 CO lines (3o3 — > 2 02 ) and (2 02 — > loi), detected with high 
signal-to-noise ratios, show a similar deviation from a Gaussian profile. They are skewed, 
with a blue-shifted peak and a weaker red line wing. The peak at v ~ 100 kms -1 may 
originate from a gas component that has been identified previously i n interferometric maps 



and may constitute the eastern part of the molecular super-shell (IWeiB et al.lll999l ). To 
better account for these line profiles, we fitted two velocity components, a "narrow" one at 
~ 100 kms -1 and a "wide" one, to each identified spectral line using the velocity and line 
width derived from the H 2 CO(3o3 — > 2 02 ) line (Table [2J Fig. H]). While the line profiles of 
the two strongest lines are fitted reliably by those two components, the result for the much 
weaker H 2 CO(3 2 i — > 2 20 ) line may be dominated by noise. 

We also inspected the spectra obtained with the other pixels of the HERA array in 
order to search for possible outflows. Filaments of SiO emission, another molecule related 
to dust-chemistry and shocks, have been found at a distance as large as 24" from the central 
plane (Fig. [2]). No significant H 2 CO emission was found in any of these spectra above a noise 
level of ~ 2mK (on the T m b scale). Note, however, that H 2 CO filaments may be present in 
the halo at locations not covered by our observations. 

In both 146 GHz spectra, the H 2 CO(2 02 — > loi) line at 145.60 GHz is partially blended 
with the HC3N(16 — > 15) line at 145.56 GHz, but the two lines can easily be separated given 
the constraints derived from the 218 GHz spectra (Table [2]). We did not detect any H 2 CO 
emission at 364 GHz from the SW lobe. The values given in Table [2] are 3a upper limits, 
derived from the rms noise with the assumption that these lines have the same width and 
central velocity as those at 218 and 146 GHz. 
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DISCUSSION 



4.1. The Physical Conditions of the Dense Molecular Gas 

For the analysis of the derived integrated line intensities, we have developed a non-LTE 
code for para-H 2 CO adopting the Large Velocity Gradient (LVG) approximation and choos- 
ing a spherically symmetric cloud geometry. The choice of a particular cloud geometry can 
affect the resulting gas densities. In the LVG approximation, the escape probability for 
a photon of a transition ij is /3y = (1 — exp(— Ty)) r^ 1 for a model with spherically sym- 
metric geometry compared to (3^ = (1 — exp(— 37y)) (Sr^) -1 for a plane-parallel geometry. 
Thus, applying a plane-parallel instead of a spherical cloud geometry, photon trapping sets 
in at lower optical depths. In the optically thick case, the particle densities d erived from 



the o bserved line ratios can be lower by up to half an ord er of m a gnitu de (e.g. IWard et al. 



20031 ). Our code is based on the collision rates with He b y I Green (1991) , the collision rates 



with H2 being approximated by a scaling factor of 1.37 (jSchoier et al 



20051 ). and includes 



41 para-H 2 CO levels, up to 210 cm -1 (300 K) above the ground state. We checked our code 
successfully by using the physical conditions derived in this study as input to the on-line 
version of RADEX0, a program using a similar non-LTE code, and by comparing its output 
with our observed and derived line ratios and optical depths. 

The velocity profiles of the two strong lines H 2 CO(2 2 — > loi) and H 2 CO(3o3 — > 2 02 ) are 
very similar in each lobe, suggesting that the H 2 CO emission is likely confined to a region 
within the smallest beam of the observations (11")- Interferometric maps of the integrated 
inte nsities of molecular transitions like 00(2 -> 1), HCO+(l -> 0) a nd HCO(F = 2 1) 
(e.g. lSeaquist et al.lll998l : IWeifi et al.ll2001bl ; lGarcia-Burillo et~aDl2002h indicate that the NE 



and the SW lobes can be approximated by circular Gaussian distributions of 7" to 8" in 
diameter (FWHM). Here we adopt a Gaussian source distribution with a size of 9 S = 7"5, 
unless noted otherwise, and convolve all intensities to a common resolution of 17". We 
calculated the expected line ratios for a large number of physical conditions, covering the 
parameter space for a kinetic temperature of Ty, n = 5 to 300 K in steps of 5 K, a molecular 
gas density of lognn2 = 3.0 to 6.0 (in cm -3 ) in steps of 0.1 and a para-H 2 CO column density 
per velocity interval of log -/V P H2co /Av = 10.5 to 14.5 (in cm -2 km -1 s) in steps of 0.1. The 
assumed ambient radiation field is the cosmic background radiation field (T bg = 2.73 K). 
Since the H 2 CO(3 22 — ► 2 2 i) line in the NE lobe is heavily blended with a nearby methanol 
line (see Section [5]) and does not add significant constraints to the physical conditions, we 
do not include the line in our LVG analysis for either lobe. 



3 http: / /www. strw.leidenuniv.nl/~moldata/ 
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Comparing the computed line ratios with our observational result (Table [3]), we can 
significantly constrain the range of possible physical conditions in the molecular lobes of 
Note that in reality, the molecular gas is likely to be inhomogeneous in its properties 
such as temperature and density. However, given the limited number of measurements, we 
are constrained to consider only a one-component model, yielding the average properties of 
the molecular gas phase traced by H 2 CO. 

The observed integrated line intensity of extragalactic sources is usually only a small 
fraction of the intensity predicted by the LVG calculations. This fraction F = T obs /T lvg , the 
so-called filling factor or dilution factor, is the product of several components F = ft, f a f v 
and provides insights into the structure of the observed source. The beam filling factor 
ft, = 6 2 s (9g + Of)^ 1 accounts for the fact that the source does not extend over the whole 
beam area. Assuming a source size of 9 X = 7'.'5, ft, = 0-163 for a beam width of Of, = 17". 
At the distance of M82 of D = 3.9 Mpc (ISakai fc Madordll999l ). a 17" beam covers an area 



with a diameter of 320 pc. Thus, it is reasonable to assume that the observed emission 
does not originate from a single molecular cloud, but rather represents an ensemble of giant 
molecular clouds, a giant molecular association. The individual clouds are expected to have 
(much) smaller radii than the overall source size and (much) smaller line widths than the 
observed beam-averaged width. Measures of this small-scale structure are the area filling 
factor f a = C r 2 R~ 2 and the velocity filling factor f v = Av c \ ou d Av# ne , where r and At> c i ou d 
are the radius and velocity width of an individual cloud, R s is the radius of the source, Avn ne 
is the observed line width, and C is the total number of clouds in the beam. Since the latter 
two filling factors can not be separated, we combine them into the small-scale dilution factor 
F sc = f a fv Note that F sc represents the volume filling factor of the velocity cube and thus 
can not be larger than unity. 

Previous studies of H 2 CO emi ssion in starburst galaxies suggest total H 2 CO abundances 



of the order of ~ 10 J to ~ 10" 8 flHiittemeister et all 119971 : IWang et all 12004 Martin et al. 



2006bl ). In the high-temperature limit, the ortho-to-p ara-H^CO ratio is olp = 3. However, 



in starburst galaxies, o/p may be closer to 1 or 2 (e.g. 
a para-H 2 CO abundance per velocity gradient of 



Hiittemeister et al.lll997h . Adopting 



A 



^pH2CO 

grad(w ) 



pH2CO 



/Av 



1 x 10 km spc, 



4 Splitting the observed H2CO line profiles in the SW lobe into two velocity components (Table [2]), the 
H2CO(2o2 — * loi)/H2CO(3o3 — > 202) line ratio suggests that the narrow component has a higher density than 
the wide component at a given abundance per velocity gradient and kinetic temperature. However, since the 
other line ratios are highly uncertain because of the low signal-to-noise ratio in the fi2CO(32i — » 220) line, 
we do not treat the two velocity component separately in the following analysis. 
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we derive the physical conditions of the molecular gas in the lobes by comparing the observed 
line ratios (Table [3l 9 S = 7".5) with those calculated in the model parameter space. While 
inter-ladder line ratios of formaldehyde are generally good tracers of the kinetic temperature 
and intra-ladder ratios sensitively trace the density at a given temperature (Section ll.il) . we 
find that at moderate gas densities of ~ 10 4 cm -3 and high kinetic temperatures of > 100 K, 
the individual line ratios are not completely independent of the other model parameters and 
that a combination of the line ratios constrains the physical properties of the molecular gas 
even more tightly. Figure shows the cuts through the parameter space where the different 
derived line ratios intersect. Also plotted are the uncertainties given in Table [3j We derive 
similar physical conditions for the two lobes, in particular a high kinetic temperature of 
~ 200 K and a moderate gas density of nH2 ~ 7.4 x 10 3 cm -3 (Table H]). Under these 
conditions, the 146 GHz line is optically thick, while the optical depth of the 218 GHz lines 
is of the order of unity and the emission at 364 GHz is optically thin. 

If the line width f c i ou d of an individual cloud is dominated by its velocity gradient 
(Aw c i ou d ~ 2rgrad(f)) and the clouds are close to virial equilibrium, a cloud gas den- 
sity of nH2 = 7.4 x 10 3 cm -3 suggests a velocity gradient of grad(f) ~ 1.0 km s -1 pc _1 , 
which is in very g ood agreement wi th the velocity gradient deduced from high-resolution 



CO observations (IWeiB et al.l l2001bl ). In this case, the para-formaldehyde abundance is 
ApH2CO ~ 1 x 10 -9 and the total formaldehyde abundance Ah2co ~ 2 to 4 x 10~ 9 , depend- 
ing on the ortho-to-para formaldehyde ratio. The mass of the molecular clouds within the 
beam of the telescope is 

M ™i = g g 2 g x io~ 19 ^ H2 ' beam ( — \ 2 ( ^ beam ^ 2 
M Q cm" 2 yMpc J \arcsec / ' 

where A^ H 2,beam = (-WpH2CO / Av) -XpH2CO At>iine F is the beam-averaged H 2 column density. 
This equation includes a factor 1.6 to accoun t for dust and other mo lecular and atomic species 



like helium. At a distance of D = 3.9 Mpc ( ISakai fc Madorelll999l ). M moi = 1.7 x 1O 8 M in 
the NE lobe and M mol = 1.4 x 10 8 M in the SW lobe. 

The small-scale dilution factor F sc = F can be derived from our LVG calculations. 
Expressed in terms of cloud properties, it is 

p _ , f _ C r 2 Av cloud 

"sc Jajv no a 

Rj Aw line 

Assuming again that the line width of an individual cloud is dominated by its velocity 
gradient, 

^-,3 R B Awi; ne F sc 

— 77 — \ — ' 

2 grad(f ) 
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yielding r = 34C~ 1 / 3 pc in the NE lobe and r = 33C~ 1//3 pc in the SW lobe for a velocity 
gradient of lkms _1 pc _1 . For comparison, the maximum radius of a cloud is given by the 
observed line width and the velocity gradient r max = At> line /(2 grad(t>)) = 66 pc and 56 pc in 
the NE and SW lobes, respectively. The fraction of the volume that the molecular clouds 
occupy within the source is 

^_ C(4/3)7rr 3 _ 2Av line F sc 
tc R 2 s l 3 grad(f ) I ' 

where I ~ 350 pc is the adopted line of sight extent of the emitting region. For a gradient of 
lkms^pc -1 , $ ~ (0.30) 3 or about 30% in each dimension of the data cube in both lobes. 
Thus, the observed volume could contain a small number of large clouds, whose diameters 
are comparable to the source size, or a large number of small clouds, but the presence of 
only a single cloud within the beam is ruled out. 



The results for other choices of A can be found in Tables [5] and [6J In the optically thin 
limit (see e.g. models Nl and N2), the small-scale dilution factor reaches its maximum value 
of F sc = 1 and the temperature-density plots are almost independent of the assumed H2CO 
column density per velocity interval. Thus, there is a lower limit to the kinetic temperatures 
and an upper limit to the gas densities consistent with our observations. We do not find a 
solution (corresponding to the intersection of the curves for the line ratios) in the calculated 
parameter space for T kin < 140 K (NE lobe) and T kin < 155 K (SW lobe), respectively. The 
upper limit to the gas density is ~ 2.5 x 10 4 in both lobes. 



The assumed source size 9 S also has a significant influence on the derived parameters 
(Table HJ). A larger source size results in a lower kinetic temperature, a higher H 2 density 
and a smaller small-scale dilution factor. But even if the H2CO emitting region is as large 
as 10", the kinetic temperat ure in the lobes is ~ 170 K, which is much higher than the dust 
temperature of T dust = 48 K ([Colbert et al.lll999l ) or the commonly assumed kinetic temper- 
atures of 30 K to 100 K. To summarize, the dominant component of the molecular gas in the 
lobes of M 82 seems to be in a warm phase of moderate density. 



4.2. Comparison with Other Excitation Studies 



Numerous previous studies have used molecular lines to derive the physical properties 
of the molecular gas in M82. Frequently, the small number of lines available to constrain 
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the free parameters has led to the ad option of a kinetic t emperature of T^- m ~ 50 K, similar 
to the dust temperature T dust = 48 K ( Colbert et al. 1999h . This assump tion is supported by 



the results of early analyz es of the low- J CO lines (e.g. 
by the ammonia study of IWeiB et al.l (j2001al . see below) 



Wild et al.lll992l ) and, more recently, 



The most comprehensive dat a base exists for CO and 13 CO lines , which have bee n 
searched for up to J = 7 -> 6 jMao et al.l boooh and J = 6 — » 5 (I Ward et al.ll2003h . 



respectively. Based on this large collection of CO lines, iMao et al.l (120001 ) confirm that a 
one-component LVG model can not fit all the observed line intensities simultaneously. Thus, 
in order to derive the properties of the high-excitation component, they restrict their LVG 
analysis to the high- J excitation lines in the submm-range, i.e. CO (7 — > 6), CO(4 — > 3) , 
CO(3 -> 2), and 13 CO(3 -> 2). In the most recent multi-line CO study, IWard et al.l (bo03h 
analyze the CO and 13 CO emission at all available transitions using a two-component plane- 
parallel LVG model. The physical properties of the high-excitation molecular gas phase 
derived from these comprehensive excitation studies are compared to our results in Table [71 
Considering that our investigation is completely independent of the other studies, the results 
are remarkably similar. The physical properties derived from our H2CO data and our LVG 
model adopting a spherically symmetric cloud morphology lie well within the ra nge of values 



found for the high-excitation component of the two-component LVG model by IWard et al. 
(120031 ). In fact, our kinetic temperature is only 20 K to 40 K higher than their median 



value, which is influenced by the limited temperature range considered, and their median 
density is lower by only a factor of ~ 2 if the effect of using a plane-parallel model versus a 
spherically symmetric one is take n into account. T he most striking difference between our 
results and the values derived by IMao et al.l (120001 ) is their lower kinetic temperature. A 
possible explanation for this discrepancy is a small contribution of the cooler low-excitation 
component to the submm CO line emission that would lower the average kinetic temperature 
of the one-component model. 

The presence of a warm gas component in nearby starburst galaxi es with T^n = 50 to 



440 K was also inferred from the detection of highly excited NH3 lines ( IMauersberger et al. 



20031 ). In this study, M82 was the notable exception in not showing evidence for a warm 
molecular gas component in its NH 3 emission. The ammonia inversion lines detected to- 
wards M82 suggest a rotational temper ature of only 29 K and thus a ki netic temperature 
of onl y ~ 60 K in the SW lobe of M82 JWeifi et alJboOlah . In contrast, iRigopoulou et al. 
(120021 ) find a warm molecular gas component with Ty, n ~ 150 K in a sample of starburst 



and Seyfert galaxies including M 82 by analyzing IR rotational H 2 emission lines. Possible 
mechanisms that could heat a large fraction of the molecular gas to tempe ratures of ~ 150 K 
include shocks, strong U V and/or X-ray irradiation (PDR, XDR) (e.g. IRigopoulou et al. 



20021 ; iFuente et al.l 120051 ) as well as the more uniform cosmic ray heating ([Bradford et al. 
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20031 ). in short processes that are likely to be found in an active environment such as a 
starburst. 

A comparison of the H 2 CO line widths and the derived molecular gas mass with the 
values derived from the CO studies suggests that the H 2 CO lines trace approximately the 
same gas as traced by the ubiquitous CO lines (Table [7]). According to our data, each lobe 
contains warm molecular gas of a bout M mo \ ~ 1.5 x 10 s Mq, which is in good agreement with 
the estimate by iMao et al.l (120001 ) of M m ? \ ~ 1 — 10 x 10 8 M w and the mass of the molecular 
ring of M mo i ~ 2.0 x 10 s M® derived by I Ward et al.l (120031 1 from the integrated intensity of 
their CO (6 — > 5) map. Having modeled the physical conditions and the conversion factor 
Xqo at 18 positions alon g the circumnuclear ring as traced by high-resolution low-J CO 



maps, 



WeiB et al.l (l2001bl ) find a molecular gas mass of M mo \ ~ 2.5 x 10 s Mq in the ring, 
whereas applying the standard conversion factor Xco = 1-6 x 10 20 cm -2 (K km s -1 ) -1 to their 
data leads to a total mass of M mo \ ~ 7.1 x 10 s Mq. 

In contrast to the results of the CO and H 2 CO investigations, the observed ammo- 
nia inversion lines seem to arise from cold, probably well shielded regions within the ISM 
( jMauersberger et al.l 120031 1 . The physical properties derived in our analysis and the CO 
studies suggest that the majority of the gas is in a warm state of comparatively low density 
(nn2 ~ 7 x 10 3 cm -3 )), which fits well with the idea that M82 is in an advanced phase of 
a starburst. In such an environ ment, ammonia is ea sily photodissociated outside of the few 



remaining well-shielded regions dFuente et al. 



NH 3 column densities reported by IWeiB et al 



993), which may explain the extremely low 



(teOOlah . 



5. METHANOL EMISSION FROM THE MOLECULAR LOBES OF M 82 

The detection of the H 2 CO(32i — > 2 2 o) line in both molecular lobes strongly sug- 
gests that th e re is H 2 CO(3 ?2 — > 2 21 ) emission at these positions too. On the other hand, 



Martin et al.l (l2006al ) very recently reported the detection of methanol emission from the 
lobes of M82 including the highly excited 5^ — ► 4^ transition. Thus, it is reasonable to 
expect CH30H(4 2 — > 3i E) emission to originate from the two lobes as well. As a re- 
sult, the observed spectral features at 218.5 GHz in both lobes most likely consist of the 
strongly blended emission of the H 2 CO(3 22 — > 2 21 ) and the CH 3 OH(4 2 — > 3iE) lines at 
218.48 GHz and at 218.44 GHz, respectively. Because of this blending, the spectral feature 
at ~ 218.5 GHz was ignored in the LVG analysis in Section H~T1 Now we can use the results 
of our LVG models to derive the expected strengths of the H 2 CO(3 22 — ► 2 2 i) emission for 
different physical conditions (Tables [5] and [6]). Interestingly, the computed values for each 
lobe vary by no more than 5% throughout the range of model parameters considered, which 
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indicates that the H2CO(322 —> 2 2 i) emission does not add significant constraints to the LVG 
analysis, probably because the similar H 2 CO(3 2 i — > 2 20 ) line is already included. In order to 
deconvolve the highly blended lines, we fixed the integrated intensity of the H 2 CO(3 22 — ► 2 21 ) 
emission to the value derived from the models and fitted the five lines with the usual con- 
straints (Table [2]). Fig. [6] shows the resulting Gaussian fits and the residual spectra in both 
lobes. The derived integrated intensities are summarized in Table [HJ 

In the SW lobe, the H 2 CO(3 22 — > 2 21 ) emission predicted by the LVG analysis agrees 
within 3cr with the previously derived integrated intensity, and the estimated contribution 
of the methanol line is well below the uncertainty of the mean of the fitted H 2 CO line. In 
the NE lobe, the methanol line appears to be about six times stronger, while the predicted 
H 2 CO(3 22 — > 2 21 ) emission is comparable to that of the SW lobe. The noise level in the 
NE residual spectrum is slightly increased at the position of the subtracted blended lines. 
However, this is also the case at the position of the subtracted H 2 CO(3o3 — > 2 02 ) line and 
may be due to a deviation of the line profile from a Gaussian shape. The result that the 
CH30H(4 2 — > 3\E) emission of the NE lobe is m uch stronger than th at of the SW lobe 



agrees very well with the line strengths reported by Martm et al.l ( 2006al ). Their intensities 



of the blended methanol lines at 97 GHz and 157 GHz are similar in both lobes. However, 
the CH 3 OH(5fe — > 4&) transitions at 242 GHz, which require high excitation levels similar to 
those needed to observe the CH30H(4 2 — > 3i E) transition, are stronger in the NE lobe than 
in the SW lobe, which might indicates a difference in excitation. 



6. CONCLUSIONS 

1. We have detected several para-H 2 CO emission lines at 146 GHz and 218 GHz in both 
molecular lobes near the center of the nearby starburst galaxy M 82. To our knowledge, 
our study constitutes the first dedicated search for H 2 CO transitions with K a > 1 and 
the first detection of H 2 CO(3 22 — > 2 21 ) emission outside of the Galactic neighborhood. 

2. Our line ratio analysis using an LVG model with spherical cloud symmetry suggests 
the presence of warm (Tkin ~ 200 K), moderately dense (n^ ~ 7 x 10 3 cm -3 ) molec- 
ular gas in the lobes near the center of the starburst activity. The ratio of column 
density to volume density and the small small-scale filling factor may indicate that the 
observed molecular gas forms large complexes of comparatively low-density molecular 
gas, possibly extended envelopes around a few remaining well-shielded cores. 

3. The physical properties of the molecular gas derived from our H 2 CO data are in very 
good agreement with parameters of the high-excitation component in the lobes of M 82 
found in recent comprehensive CO studies, but differ from the results of a multi-line 
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NH3 investigation. Thus, the para-h^CO lines seem to trace the probably dominant 
high-excitation component of the molecular gas in M82 very well, while the NH 3 
emission might originate predominantly from cold cloud cores. 

4. The total mass of the molecular gas observed in the lobes is a few times 10 s M Q , similar 
to the mass determinations of other molecular line studies. 

5. The selected para-H^CO transitions are good tracers of the molecular gas in starburst 
galaxies, even if a large fraction of the gas is in a warm phase (Tkm ~ 150 K) of 
moderate density (n m < 10 4 cm~ 3 ). 

6. We find strong evidence for CH 3 OH(4 2 — > 3i E) emission in the spectrum of the NE 
lobe at 218 GHz and also estimate the intensity of this e mission line in the SW lobe. 



We thus confirm the detection of methanol in M 82 by iMartm et all (l2006bl ) with 



measurements of another transition line that are consistent with the earlier findings. 
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Table 1. Summary of the observed para- (p) and ortho-H 2 CO (o) transitions. 



Transition 




i//[GHz] 


e b */["] 


rms b / [mK] 


H 2 CO(2 02 


- loi) 


P 


145.6029 


17 


1.9, 1.8 


H 2 CO(3 03 


^ 2 02 ) 


P 


218.2222 


11 


2.7, 2.4 


H2CO(322 


->2 2 i) 


P 


218.4756 


11 


2.7, 2.4 


H 2 CO(3 2 i 


^ 2 20 ) 


P 


218.7601 


11 


2.7, 2.4 


H 2 CO(5 24 


4 23 ) 


P 


363.9459 


13 


5.4, n/a 


H 2 CO(5 4 2 


-44l) 


P 


364.1032 


13 


5.4, n/a 


H 2 CO(5 4 i 


^ 4 40 ) 


P 


364.1032 


13 


5.4, n/a 


H 2 CO(5 33 


4 32 ) 


o 


364.2751 


13 


5.4, n/a 


H 2 CO(5 32 


^4 31 ) 


o 


364.2889 


13 


5.4, n/a 



a half power beam width (FWHM) 

b rms in T m b at lOkms -1 resolution at the positions of the 
SW lobe and the NE lobe, respectively 
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Table 2. Input parameters and results of the Gaussian fits to the observed emission lines. 



Parameter 




NE lobe 


SW lobe 


SW lobe/narrow 


SW lobe/ wide 


vq [kms -1 ] 




301.0 ±3.6 


132.2 ±4.8 


106.3 ±5.8 


174.3 ± 12.7 


Awiino [kms -1 ] 




130.7 ±4.1 


111.2 ±6.4 


48.3 ± 10.5 


77.1 ± 18.6 


/ T mb <fo(H 2 CO(2 02 - 


-loi)) 


3.69 ±0.10 


2.76 ±0.10 


1.34 ±0.07 


1.60 ±0.15 


/T mb ^(HC 3 N(16- 


•15)) 


0.46 ±0.11 


0.70 ±0.09 


< 0.19 


0.42 ±0.08 


/T mb dw(H 2 CO(3o3- 


- 202)) 


2.75 ±0.18 (2.73) 


2.09 ±0.16 (2.09) 


0.98 ±0.08 


1.10 ±0.12 


/T mb( i«(H 2 CO(3 22 - 


-2 21 )) 


0.05±0.17 (0.11) 


1.09±0.15 (1.10) 


0.57±0.08 


o.4i±o.io 


/ T mbC fe(H 2 CO(3 21 - 


- 2ao)) 


0.67 ±0.15 (0.68) 


0.53 ±0.14 (0.53) 


0.25 ±0.09 


0.21 ±0.08 


/ T mb cfe(CH 3 OH(4 2 




1.45±0.21 (1.39) 


(< 0.27) 






Jr mb d W (HC 3 N(24- 


•23)) 


(0.17) 


(< 0.27) 






/T mb dv(H 2 CO(524 - 


- 4 23 )) 




< 0.54 






fT mh dv(R 2 CO{5 4x - 


- 4 4 ,)) a 




< 0.54 






jT mb dv(U 2 CO(5 3x - 


- 4 3y )) b 




< 0.54 







integrated intensity of the highly blended lines H 2 CO(54 2 — > 441) and H 2 CO(54i — ► 440) 
integrated intensity of the highly blended lines H 2 CO(5 33 — > 4 32 ) and H 2 CO(5 32 — > 4 3 i) 



Note. — All integrated intensities in Kkms" 1 (main-beam temperature scale) and before correction for beam 
width. The last two columns show the line intensities in the SW lobe that result from the assumption of two 
velocity components ("narrow" and "wide"). At each position, the velocity vo, the line width At>ii ne and their 
errors were derived by simultaneously fitting Gaussians with a fixed frequency offset and the same (unknown) 
line width to the detected lines in the 218 GHz spectrum (NE lobe: 4 components, SW lobe: 3 components, SW 
lobe/narrow and wide: 2 components to the H 2 CO(3o 3 — » 2 02 ) line alone). The derived velocities and line widths 
were then fixed for the fits to the other spectra. Values in parentheses show the results of a simultaneous fit to all 
five potentially present lines in the 218 GHz spectra (H 2 CO(3 03 -» 2 02 ), HC 3 N(24 -> 23), CH 3 OH(4 2 -> 3i E), 
H 2 CO(3 22 — > 2 2 i), and H 2 CO(3 2 i — > 2 20 )); values in italics indicate highly blended lines, whose fit is uncertain; 
upper limits in this table are 3a limits to the integrated intensities of the spectral lines using the adopted line 
widths. 
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Table 3. Beam-corrected line ratios derived from the observed line intensities and the 

assumed source size 6 S 



line intensity ratio / source, 6> s a NE, 5" NE, 775 NE, 10" SW, 5" SW, 775 SW, 10' 



H 2 CO(2 02 - loi)/H 2 CO(3 3 - 2 02 ) 
H 2 CO(2 02 - l i)/H 2 CO(3 21 -» 2 20 ) 
H 2 CO(3 03 -» 2 02 )/H 2 CO(3 2 i -» 2 20 ) 
H 2 CO(2 02 -» l i)/H 2 CO(5 24 -» 4 23 ) 



2i 



3+0.20 
3 -0.22 

11.91±1;?1 

4.i3±? : f§ 



2.6il8:5§ 

10.79™ 



4.13 



+0.80 
1.12 



0.16 
0.18 
+ 1.80 
-2.54 
4 1Q+0.80 



2.36 
9.75 



2.84lg;» 
11. 26±|" li 

3.96+?; 8 2i 
> 8.28 



9 KS+0.21 
Z.OO_ 24 

10 20+ 2 09 

1U.ZU_ 3 QQ 

3.96±?:I5 

> 7.84 



O OO + 0.19 

Z.OO_ 2i 

9.22i 2 ' 8 7 ? 
3.96t?;l? 
> 7.40 



a selected lobe, assumed source size 9, 
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Table 4. Results of our LVG calculations and the derived physical properties assuming 





A 

A = 


1 x 10 y km 


spc. 








parameter / source, 9 s a 


NE, 5" 


NE, 7'.'5 


NE, 10" 


SW, 5" 


SW, 7'.'5 


SW, 10" 


T kin [K] 


230 


191 


161 


256 


209 


181 


logriH2 [cm -3 ] 


3.67 


3.87 


4.04 


3.67 


3.87 


4.02 


log^V p H2Co/Aw b [cm -2 km" 1 s] 


13.16 


13.36 


13.53 


13.16 


13.36 


13.51 


7145.6° 


2.508 


3.628 


4.306 


2.473 


3.551 


4.173 


T218.2° 


0.636 


1.345 


1.933 


0.642 


1.355 


1.943 


T 218.5° 


0.261 


0.432 


0.520 


0.274 


0.455 


0.557 


T218.8° 


0.153 


0.273 


0.343 


0.163 


0.291 


0.368 


^364.0° 


0.003 


0.007 


0.010 


0.004 


0.008 


0.011 


pc,d 


0.025 


0.019 


0.013 


0.037 


0.018 


0.013 


P c,e 
sc 


0.24 


0.11 


0.08 


0.23 


0.11 


0.08 


P = n H2 T ki n [10 6 cm- 3 K] 


1.1 


1.4 


1.8 


1.2 


1.6 


1.9 


A r pH2CO,bcam f [10 13 Cm" 2 ] 


4.7 


5.7 


5.8 


6.0 


4.6 


4.7 


A r H2,bcam S [10 22 Cm" 2 ] 


4.7 


5.7 


5.8 


6.0 


4.6 


4.7 


M mol h [1O 8 M ] 


1.4 


1.7 


1.7 


1.8 


1.4 


1.4 



a selected lobe, assumed source size 9 S 
b calculated from A and logriH 2 
c evaluated at the nearest grid point 

d dilution factor F = T b s /T\ vg in a 17" beam; the 218 GHz data are scaled to a 17" beam assuming 
a source size of 9 S 

c small-scale dilution factor F sc = F/ 6 _1 = f a f v = F/0.163 for 9 S =7'/5. 
f beam-averaged column density iVp H 2CO,bcam = (^ P H2Co/A«) Avu ne F 
s beam-averaged column density, assuming X p h 2 co = 1 x 10~ 9 
h molecular gas mass, assuming X p H2CO = 1 x 10~ 9 and D = 3.9 Mpc 



Table 5. Results of our LVG calculations for various physical conditions of the molecular gas in the NE lobe. 



model 


-Man 


log 7JH2 


log iV n H2ro / Au 


X n H?m / grad v 


Ti46 a 


T218.2 a 


T 218.5 a 


T"218.8 a 


r 36 4 a 


pa.b 


F sc a ' c 


f Tmh dv^ia c; a ' 
j mu zio.o 




[Kl 




[cm - 2 km - 1 s] 


[lO^km" 1 spc] 
















[Kkms" 1 ] 


Tl 


150 


4.24 


12.73 


0.100 


0.78 


0.22 


0.09 


0.06 


0.002 


0.059 


0.36 


0.53 


T2 


175 


4.02 


13.18 


0.469 


2.41 


0.79 


0.29 


0.18 


0.006 


0.027 


0.16 


0.54 


T3 


200 


3.79 


13.47 


1.55 


4.72 


1.81 


0.53 


0.33 


0.008 


0.020 


0.12 


0.55 


T4 


225 


3.64 


13.61 


3.03 


6.24 


2.23 


0.58 


0.34 


0.007 


0.023 


0.14 


0.55 


T5 


250 


3.52 


13.75 


5.51 


9.52 


4.02 


0.84 


0.51 


0.010 


0.017 


0.11 


0.54 


T6 


275 


3.40 


13.88 


9.81 


11.96 


5.24 


0.95 


0.59 


0.010 


0.016 


0.10 


0.55 


T7 


300 


3.28 


13.99 


16.65 


15.07 


6.82 


1.11 


0.67 


0.011 


0.016 


0.10 


0.55 


Nl 


142 


4.38 


12.00 


0.0135 


0.15 


0.04 


0.19 


0.13 


0.0007 


0.202 


1.24 


n/a e 


N2 


146 


4.31 


12.50 


0.0503 


0.49 


0.14 


0.06 


0.04 


0.002 


0.078 


0.48 


0.52 


N3 


159 


4.13 


13.00 


0.2407 


1.54 


0.46 


0.19 


0.12 


0.004 


0.036 


0.22 


0.53 


N4 


206 


3.76 


13.50 


1.7842 


4.68 


1.81 


0.53 


0.33 


0.008 


0.020 


0.12 


0.55 


N5 


298 


3.28 


14.00 


17.0392 


15.04 


6.80 


1.11 


0.67 


0.011 


0.016 


0.10 


0.55 



Note. — For the T models, a kinetic temperature was specified, while the N models used a specific para-H 2 CO column density per velocity 
interval as an input parameter. 

Evaluated at the nearest grid point 

b dilution factor F = T bs/7ivg in a 17" beam; the 218 GHz data are scaled to a 17" beam assuming a source size of 8 S =7"5 
c small-scale dilution factor F sc = F/ 6 _1 = f a f v = F/0.163 for 9 S =7'/5. 
Calculated integrated intensity of the H 2 CO(32i — > 2 2 o) line 
e value ambiguous 



Table 6. Results of our LVG calculations for various physical conditions of the molecular gas in the SW lobe. 



model 




log 7lH9 


log N„u9m 1 Aw 


XnH9rn / £rad v 


7"l46 a 


7"218 2 a 


Tola t; a 


T91S S a 


T-?fi4 a 


pa,b 


p a,c 
sc 


J IllD " /t/ Zl?S.O 




[K] 




[cm -2 km - 1 s] 


[lO^km" 1 spc] 
















[Kkms^ 1 ] 


T2 


175 


4.17 


12.83 


0.148 


0.943 


0.287 


0.130 


0.085 


0.004 


0.038 


0.23 


0.42 


T3 


200 


3.93 


13.28 


0.727 


3.014 


1.024 


0.366 


0.230 


0.006 


0.022 


0.13 


0.43 


T4 


225 


3.78 


13.45 


1.52 


4.605 


1.832 


0.561 


0.354 


0.009 


0.016 


0.10 


0.43 


T5 


250 


3.65 


13.60 


2.89 


5.790 


2.393 


0.671 


0.423 


0.010 


0.015 


0.09 


0.43 


T6 


275 


3.53 


13.72 


5.03 


7.726 


2.948 


0.717 


0.432 


0.008 


0.018 


0.11 


0.43 


T7 


300 


3.43 


13.86 


8.74 


11.78 


5.282 


1.022 


0.625 


0.011 


0.014 


0.08 


0.43 


Nl 


159 


4.34 


12.00 


0.0148 


0.153 


0.041 


0.020 


0.014 


0.0006 


0.203 


1.25 


0.41 


N2 


164 


4.27 


12.50 


0.0551 


0.472 


0.140 


0.065 


0.044 


0.002 


0.064 


0.39 


0.42 


N3 


183 


4.10 


13.00 


0.2579 


1.492 


0.470 


0.201 


0.130 


0.005 


0.029 


0.18 


0.43 


N4 


230 


3.75 


13.50 


1.8258 


4.575 


1.827 


0.564 


0.358 


0.009 


0.016 


0.10 


0.43 



Note. — For the T models, a kinetic temperature was specified, while the N models used a specific para-H 2 CO column density per velocity 
interval as an input parameter. 

Evaluated at the nearest grid point 

b dilution factor F = Tobs/Tivg in a 17" beam; the 218 GHz data are scaled to a 17" beam assuming a source size 9 S =7"5 

c small-scale dilution factor F sc = F f^ 1 = f a f v 

Calculated integrated intensity of the H 2 CO(32i — » 2 2 o) line 



Table 7. Comparison of our results to those of other excitation studies. 



parameter /study 


our data (NE/SW) 


Ward/median (NE/SW) a - b 


Ward/range (NE/SW) a ' c 


Mao (NE/SW) d 


9b [arcsec] 


17 


24.4 


24.4 


22 


Awii no [kms^ 1 ] 


131/111 


180/160 


180/160 


- 120/- 130° 


T kin [K] 


191/209 


170/170 


> 50/> 50 


60 - 130 


nm [10 3 cm~ 3 ] 


7.4/7.4 


0.6/1.0 


0.3-100/0.3-20 


2.0-7.9 


A r H2,bcam [10 22 Cm- 2 ] 


5.7/4.6 


1.3/1.0 


0.08-5.0/0.16-12.6 


~ 10 


M mol [10 8 M ] 


1.7/1.4 


2.0 f ' 9 




1 - 10 s 



a high-excitation component of a two-component LVG model assuming a plane-parallel geometry, using all available 
CO and 13 CO lines; for a spheri cal model, riH2 ~ 5 times larger, Nm ~ 2 times larger and Tki n better constrained and 
slightly lower (jWard et al.ll2003h 

b median likelihood of possible values, not a self-consistent solution 

c full range of possible values (95% confidence interval) 

d one-component LVG model assuming a spherically s ymmetric cloud geometry, using only the high-excitation lines 
CO(7 -> 6), CO(4 -> 3), CO(3 -> 2), and 13 CO(3 -> 2) {iVlao et al J 1 200(1 ) 

C FWHM, derived from CO (2 -> 1) and CO(l -> 0) spectra 

f using the average CO column density derived from the integrated CO (6 — > 5) intensity map 
s scaled to a distance of D = 3.9 Mpc, assumed to be the total mass 
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Table 8. Derived integrated intensities of the emission lines at 218.5 GHz at 11" 

resolution. 





NE lobe 


SW lobe 


/T mb <MH 2 CO(3 22 2 21 )) modcls [Kkms- 1 ] 


1.07 


0.84 


/T mb ^(CH 3 OH(4 2 3i £?)) fitted [Kkms- 1 ] 


0.61 


0.10 


3cr [Kkms" 1 ] 


0.26 


0.27 


residual noise level at 218.5 GHz [mK] 


3.0 a 


2.3 


rms noise level over the whole baseline [mK] 


2.4 


2.7 



a for comparison, the residual noise level at the position of the 
H 2 CO(3 03 2 02 ) line is 3.3 mK 
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Fig. 1. — Contours of the HCO(F = 2 — ► 1) emission supe rposed on a high-resolutio n 
CO (2 — ► 1) map of the molecular ring in M82 adapted from iGarcia-BuriUo et al.l (120021 ). 
Map offsets are relative to the dynamical center of the galaxy (a 2 ooo = 09 h 55 m 51!9, <5 2 ooo = 
69°40'47'.'l) with the x-axis along the plane of the molecular ring (P.A.= 70° east of north). 
The observed pointing positions and the corresponding beam widths are marked by the solid 
(11" at 218 GHz) and the dashed (17" at 146 GHz) circles. 
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Fig. 2. — Contours of the velocity-integrated SiO(t> = 0, J = 2 — 1) emission sup e rpose d 
on a radio continuum emission image at 4.8 GHz, adapted from iGarcia-Burillo et al.l (120011 ). 
Map offsets are relative to the dynamical center of the galaxy («2ooo = 09 h 55 m 51!9, <52ooo = 
69 o 40'47'.'l). The large circle delimits the primary beam field of the SiO interferometric 
observations at 87 GHz (55"), while the synthesized beam is shown in the bottom left corner. 
The white square marks the position of SNR 41.95+57.5. The radio continuum filament is 
indicated by an arrow. The pointing positions observed with the HERA array and the beam 
width of 11" are marked by the gray and black solid circles, where the gray circles cover the 
molecular lobes in M82. 
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Fig. 3. — Observed spectra of the NE (left) and the SW lobe (right) of M82. Each spectrum 
is labeled with the frequency the receiver was tuned to. Thus, the velocity scale of the 
218.48 GHz spectra refers to the H 2 CO(3 22 -> 2 21 ) line. The H 2 CO(3 03 -> 2 02 ) and the 
H 2 CO(3 2 i — > 2 20 ) transitions are offset by 348.5kms~ 1 and — 389.7km.s~ 1 , respectively, 
while the CH 3 OH(4 2 -> 3 1 E) line is offset by 49.4 km s" 1 . In the 146 GHz spectra, the 
HC 3 N(16 — > 15) line is offset by 86.46 kms -1 from the H 2 CO(2 02 — > loi) emission. Gaussian 
line profiles were fitted to every individual line assuming the same line width and velocity for 
all lines. All identified lines of a spectrum were fitted simultaneously with these constraints. 
The curves show the Gaussian fit to each individual line as well as the spectrum resulting 
from a superposition of all identified lines: a = H 2 CO(2 02 — > loi), b = HC 3 N(16 — > 15), c = 
H 2 CO(3 21 -> 2 20 ), d = H 2 CO(3 22 -> 2 21 ), e = CH 3 OH(4 2 -> 3i E), f = H 2 CO(3 03 -> 2 02 ), g 
= H 2 CO(5 33 -> 4 32 ) + H 2 CO(5 32 -> 4 31 ), h = H 2 CO(5 24 -> 4 23 ). 
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Fig. 4.— Spectra of the SW lobe observed at 146 GHz (left) and 218 GHz (right). To 
better match the observed line profiles, two velocity components (v\ and i> 2 ) with a mutual 
separation of GSkms -1 were fitted to each previously identified line. The curves show the 
Gaussian fit to each individual line as well as the spectrum resulting from a superposition of 
all identified lines: a = H 2 CO(2 02 -> l 01 ) wl , b = H 2 CO(2 02 -> l i)«2, c = HC 3 N(16 -> 15)„i, 
d = HC 3 N(16 -> 15)„ 2 , e = H 2 CO(3 21 -> 2 20 U f = H 2 CO(3 21 -> 2 20 )„ 2 , g = H 2 CO(3 22 -> 
2 2 i)„i, h = H 2 CO(3 22 -> 2 2 i)„i, i = H 2 CO(3 03 -> 2 02 )„i, j = H 2 CO(3 03 -> 2 02 ), 2 . Note 
that the H 2 CO(2 02 — > loi)v2 and the HC 3 N(16 — > 15)„i lines are strongly blended with a 
separation in velocity of only 16.5 km s -1 . 
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Fig. 5. — Results of our LVG analysis for the NE and SW lobes shown as cuts through the 
3-D parameter space (T kin in K, n H2 in cm -3 , N pm co/Av i n cm" 2 km -1 s) along a plane 
of constant para-H 2 CO column density per velocity interval (top) and of constant kinetic 
temperature (bottom). The thick lines represent the following line ratios: H2CO(3o3 — > 
2 02 )/H 2 CO(3 21 - 2 20 ) (solid), H 2 CO(2 02 - 1 01 )/H 2 CO(3 03 - 2 02 ) (dashed), H 2 CO(2 02 - 
loi)/H 2 CO(3 2 i — > 2 20 ) (dotted). The thin lines outline the corresponding uncertainties given 
in Table El The dashed-dotted line in the plots for the SW lobe denotes the area that is 
ruled out by the lower limit to the line ratio of H 2 CO(2 02 — > loi)/H 2 CO(5 2 4 — > 4 23 ). The 
adopted source size is 9 S = T!5 and the assumed para-H 2 CO abundance per velocity gradient 
is A = 1 x 1CT 9 km -1 s pc. 
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Fig. 6. — Top: Gaussian fits derived from the LVG analysis (Section [5]) superposed on 
the 218 GHz spectra of the NE lobe (left) and the SW lobe (right); the curves show the 
Gaussian fit to each individual line as well as the spectrum resulting from a superposition of 
all identified lines: a = H 2 CO(3 2 i -> 2 20 ), b = H 2 CO(3 22 -»■ 2 21 ), c = CH 3 OH(4 2 -> 3i£), 
d = H 2 CO(3q3 — > 2q 2 ). Bottom: Residual spectra after subtraction of the Gaussian fits. 



